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SOLITONS OF DISCRETE CURVE SHORTENING 


CHRISTINE RADEMACHER AND HANS-BERT RADEMACHER 


Abstract. Eor a polygon x = in R" we consider the polygon 

{T{x))j = {xj-i + 2xj + Xj+i} /4. This transformation is obtained by 
applying the midpoints polygon construction twice. Eor a closed polygon 
or a polygon with finite vertices this is a curve shortening process. We 
call a polygon x a soliton of the transformation T if the polygon T{x) 
is an affine image of x. We describe a large class of solitons for the 
transformation T by considering smooth curves c which are solutions of 
the differential equation c{t) = Bc{t) + d for a real matrix B and a vector 
d. The solutions of this differential equation can be written in terms of 
power series in the matrix B. For a solution c and for any s > 0, a € R 
the polygon x{a,s) = (xj(a, s)j) j^z; xj (a, s) = c(a + sj) is a soliton of 
T. For example we obtain solitons lying on spiral curves which under 
the transformation T rotate and shrink. 


1. Introduction 

For an infinite polygon {xj)j£i given by the vertices Xj in the vector space M"' 
the midpoints polygon is defined by M{x) where M{x)j := (xj +Xj+i)/2; j E 
Z. If the polygon is closed or rather periodic, i.e. if for some N : xj+n = xj 
for all j € Z then this midpoint mapping M defines a curve shortening pro¬ 
cess, i.e. the polygon M(x) is shorter than the polygon x unless it is a single 
point. If we iterate the process for a closed polygon {xj)j=i^,„^N it converges 
to the barycenter (xi + ... + X]sr)/N. This elementary construction was al¬ 
ready used by Darboux in 1878. He also showed in [6] that in the plane in 
the general case the sequence M^{x)/ cos(7r/V)^ converges for /c —>• oo to an 
ellipse. These results for polygons in Euclidean space were later rediscovered 
and extended by several authors, for example Kasner m, Schoeneberg m, 
and Berlekamp et al. [2]. Bruckstein & Shaked [3] discuss the relation with 
iterative smoothing procedures in shape analysis and recognition. Nowadays 
applications of discrete curve shortening are discussed in several papers. For 
example Smith et al. |15| present its connection with the rendezvous problem 
for mobile autonomous robots. 

We modify the curve shortening process M for an infinite polygon (i.e. not 
necessarily closed) in Euclidean space as follows: Instead of the midpoint 
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mapping M we apply the midpoint mapping twice and use an index shift, 
i.e. we define the polygon T{x) = {T{x)j)j£z by the Equation 

(1) T{x)j = I {xj-i + 2xj + Xj+i} , 

i.e. T{x)j = . Introducing the index shift has the following 

advantage: Since 

( 2 ) T{x)j - Xj = ^ + Xj+i} 

we can view this process T as a discrete approximation of the semidiscrete 
flow s i-A Xj{s) defined by 

dxj(s) 

(3) ^ Xj-i{s) - 2xj[s) + Xj+i(s), Xj{0) = Xj 

on the space of polygons. On the space of closed polygons this flow is the 
negative gradient flow of the functional F 2 , cf. Section |H Semidiscrete flows 
are discussed for example by Chow&: Glickenstein [5]. 

The polygon T{x) is formed by the midpoints T{x)j of the medians through 
Xj of the triangle formed by Xj-i,Xj,Xjj^i. Linear polygons Xj = ju + 
v; u,v & M"" are the hxed points of T. 

Since the process T is affinely invariant it is natural to consider polygons 
which are mapped under T onto an affine image of themselves. We also 
use the term soliton for these polygons in analogy to the case of smooth 
curves or manifolds which are mapped under the mean curvature flow onto 
the same curve or rather manifold up to an isometry, cf. Hungerbiihler 
&: Smoczyk |9], Hungerbiihler &: Roost |8] and Altschuler |T]. We call a 
polygon X = ixj)j a soliton of the process T if there exists an affine mapping 
X G M"" I—7> Ax + b £ M” for a matrix A and a vector b such that 

(4) T{x)j = ^ {xj-i + 2xj + Xj+i} = Axj + b 

for all j € Z. The main idea of this paper is to consider not only polygons, 
but also smooth curves c which are affinely invariant under an analogous 
process on curves: We adapt the process T to curves in the following sense: 
We associate the one-parameter family c*, s G M to the smooth curve c : 

(5) Cs{t) := ^ {c{t - s) + 2c{t) + c{t + s)} . 

For s > 0,a G M we define the polygon x{a,s) = {xj{a,s)j)j^^-,Xj{a,s) = 
c(a -|- sj). Then the smooth curve c defines solitons of the form x(a,s) for 
the curve shortening process T if the following holds: For some e > 0 and 
any s G (0, e) there is a one-parameter family of affine maps x G M” i-A 
A(s)x -|- b{s) G M” such that for any s G (0, e) and all t G M : 

(*) Cs{t) = ^ {c(f - s) + 2c{t) + c{t + s)} = A{s)c(t) + b{s). 

In Figure [T] we show a polygon x = (xj) which is a soliton of the pro- 
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Figure 1. polygon 
as soliton of T 


Figure 2. smooth 
curve as soliton of T 


cess T. Its vertices Xj = c(0.4-j),j G Z lie on the smooth curve c{t) = 
(cos(2t), cos(3t)), which is also a soliton of the process T, see Figure [2j 
Denote by A the diagonal matrix with entries (1 + cos( 0 . 8))/2 and (1 + 
cos(1.2))/2. Then in Figured] the polyon x and its image T[x) are shown, 
which satisfy Equation (|d|). In Figure [2| the smooth curves c and co ,4 are 
shown, which satisfy Equation Q for s = 0.4, ^(0.4) = A, 6(0.4) = 0. Hence 
the process T in this case corresponds to a scaling. This example belongs to 
Case (la) in Section [5] Here H(0.4) is given by Equation dHj) where B is the 
diagonal matrix with entries 6 i = —4 and 62 = —9. 

It is the main result of this paper that solutions of Equation (|*j), i.e. smooth 
curves c defining solitons for the curve shortening process T, can be charac¬ 
terized as solutions of the inhomogeneous linear differential equation 

( 6 ) c = Bc{t) + d 

of second order with constant coefficients with B = 2 H"(0), d = 2 6"(0), cf. 
Theorem dl and Theorem [2] The solitons c and the maps ^(s) or rather h[s) 
can be described in terms of the power series 

k —0 k —0 

cf. Proposition dl and Proposition d) in Section [2] . 

We will show, that for any real matrix B solutions of Equation Q are solitons 
of the curve shortening process and the matrix A(s) of Equation Q is given 

by 

( 8 ) H(s) = ^{ 1 -hcoij(s)} . 

The vectors b{s) depend on the structure of the matrix B and vanish for 
closed solitons. 

The one-parameter family s e->■ Cs{t) defined by Equation Q associated to 
a soliton c can be used to define an affinely invariant solution of the wave 
equation, cf. Remark [3l The relation of polygonal curve shortening and the 
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curve shortening flow for smooth curves is discussed in many papers, cf. m 
Sec.3], a reference for curve shortening flows is the book Self-similar 
solutions of the Euclidean curve shortening flow in the plane, i.e. the mean 
curvature flow for curves in Euclidean space, are discussed by Halldorsson [7], 
see also Hungerbiihler & Smoczyk j9] and Altschuler [T] . These questions lead 
to systems of non-linear ordinary differential equations. 

In Section 0] we show that the solitons of the curve shortening T are also 
solitons of the semidiscrete flow defined in Equation ([3]). In Section [5] we 
discuss the planar case n = 2 in detail and present examples. For example 
the parabola is a soliton for which the curve shortening leads to a translation. 
We also obtain various spirals as solitons which rotate and shrink and closed 
curves of Lissajous type, which scale under the mapping T. The results of 
Section [5] should be compared with the zoo of solitons obtained for the Eu¬ 
clidean curve shortening flow by Halldorsson [7], see also Hungerbiihler & 
Smoczyk |9] and Altschuler [T]. 


2. System of linear differential equations of second order 


Let B G M(re;M) be a n x n matrix with real entries. We denote by 1 the 
identity matrix and use the following notation: 

Using the matrix exponential exp(H) := ^'^oB^/kl one can also define the 
matrix functions cosh(H),cos(H),sin(H),sinh(H) as power series in B. For 
commuting matrices B,C (i.e. BC = CB) we have 

exp(i? + C) = exp(H) exp(C') = exp(C') exp(H) , 

cf. m ch.6.1]. 


Definition 1. For a matrix B G M(n,M) we define the following mappings 
COB, sis : M —)■ M (n, M) : 


cosit) 

siB{t) 


°° £2k 


E 

k=0 


{2k) 


J.2 

= 1 + —B + -B'^ + ... 


2V 


4 !' 


E 






Bk+l 

- -B^ = tl+'-B+ '-B^ + 

^^^(2A: + 1)! 3! 5! 


These power series are obviously convergent and we obtain the following 


Proposition 1 (Properties of cob{ t), siFor a matrix B the above 
defined mappings coB(t), siB(t) satisfy the following (differential) equations: 

^coB(t) = co^it) = B ■ siBit ); sifi{t) = coB{t) 

co'f{t) = B ■ coB{t ); si'f{t) = B ■ siB{t) 

COB{t)^ — B ■ siB(t)^ = 1. 


Given a real number 6 G M we denote by cos6,sinb : M —>■ M the unique 
solutions of the differential equation f" = bf with cosf,(0) = l,cos^(0) = 
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0, sint,(0) = 0,smj(0) = 1. Hence 
(9) 



1 ; 6 = 0 
cos ; 6 < 0 


6 > 0 

6 = 0 ; sinb(t) = < 


/ 


sinh /Vb ; 6 > 0 
t ] b = 0 


sin (-v/^t) /; 6 < 0 


Then we obtain for a real number 6 and the matrix B = bl : 

coB{t) = cobi{t) = cosfe(t)l; siB(t) = sibi{t) = sinb(t)l. 

If H C G M {n, C) then cosit) = cosh(Ct), C ■ sisit) = sinh(fC). For 

B = —C^, C G M {n, C) we obtain cosit) = cos(C't), C ■ siB{t) = sin(C't). 

Proposition 2 (Homogeneous differential equation). For a matrix 
B G M(n; M) and vectors v,w G M"' the linear system of ordinary differential 
equations (with constant coefficients) of second order: 

(10) c{t) = Bc{t) 

with initial values c(0) = v, c(0) = w has the unique solution 

(11) c{t) = c[v,w]{t) := coB{t){v)+ siB{t){w). 

Then the one-parameter family of curves s i-A c* for s G M defined by Equa¬ 
tion dSD satisfies 


( 12 ) 


Cs{t) = A{s) c{t) 


with 



i. e. Cs is a linear image of c for any s and Equations (|*]) and Equations ([8]) 
hold with b{s) = 0 . 

Proof. From Proposition [T] it follows that c{t) = c[v,w]{t) is a solution of 
Equation (fTOl) with c(0) = v, c(0) = w. Equation ([5]) implies that Cg is the 
unique solution of Equation (fTOll with initial conditions: 


c^(0) = i{c(-s)+ 2c(0)+ c(s)} = i{l+ cob(s)} (u) 
c^(0) = |{c(-s)+ 2c(0)+ c(s)} = ^{1+ coij(s)} (rr). 


Using Equation dill) we can write 



□ 
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Remark 1 (Addition rules). The second order linear system ffT0)l of 
ordinary differential equations is equivalent to the following first order linear 
system of ordinary differential equations 

— ( ^ ° ^ ^ \ 
dt V c(t) J \ B 0 J \ c{t) ) 

with constant coefficients. The solution of this system with initial values 
c(0) = u, c(0) = w is given by 


c(t) 

c{t) 


= exp 


0 1 
B 0 


V 

w 


COB{t) sisit) W A 

B ■ SlBit) COB{t) / V ^ / 


Hence the curve t i-A (c(t),c(t)) is the orbit of the point {v,w) G M” © M"" 
under a one-parameter group of linear transformations. 

Since 


exp 


0 = cos(t + s) siB{t + s) \ 

B 0 ’) \B-siB{t + s) coB{t + s) J 

= ®p((b 0 )*)■“'’(( s 0 

^ / cos(t) sis(t) \ _ f CO b{s) SlBis) \ 

\B-siB{t) COB{t) J \B-siB{s) COb{s) j 



we obtain the following addition rules: 


(13) COB{t + s) = COB{t) ■ cob{s) + B ■ siB{t) ■ siB{s) 

(14) siB{t + s) = COB{t) ■ SlBis) +siB{t) ■ cob{s) . 


Remark 2 (Roots of B). If there is a matrix C G M(n;M), such that 
B = then the curve c{t) = c[v,Cv]{t) = exp(tC')(u) satisfies c{t + s) = 
exp(sC)c(t), i.e. the curve is the orbit of a point under a one-parameter 
group of affine transformations in M”'. 


Let c be a solution of the inhomogeneous linear differential equation c(t) = 
Bc{t) + d, let U G M(n;M) be an invertible matrix and e a vector, then 
f{t) = Uc{t) + e solves the differential equation f{t) = UBU~^f{t) + 
(l — UBU~^^ e + Ud. Therefore one can assume that B G M(n;M) has 
already complex Jordan normal form cf. m Thm. 5.4.10] or real Jordan 
normal form, cf. im Thm. 5.6.3]. If an eigenvalue A of A is not real then the 
conjugate value A is an eigenvalue, too. In Section [5] we discuss the possible 
real Jordan normal forms for n = 2. For a complex number A the Jordan 
block JmW is given by: 




/ 

A 

1 

0 

0 . 

■ ° ^ 





0 

A 

1 

0 . 

. 0 


(15) 

Jm{^) • — 







> 




0 

0 

0 

0 . 

. 1 




V 

0 

0 

0 

0 . 

• 
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The complex Jordan normal form of a real matrix B consists of Jordan blocks 
of the form > 1 for an eigenvalue A G C. 

The real Jordan normal form of a real matrix B consists of two different 
types of Jordan blocks: For a real eigenvalue A there are Jordan blocks of 
the form Jm(A),m > 1. If /i = a + z/3,a G M,/3 G M \ {0} is a non-real 
eigenvalue, then for some m > 1 real Jordan blocks of the form 

(16) (o, /3) := ( ) } 2 m, 

occur. The real Jordan block J 2 m {oi, (3) is always invertible. Hence the 
Jordan normal form of a singular real matrix contains a nilpotent Jordan 
block of the form Jm{0),m > 1. Therefore it is sufficient to discuss the 
following cases: 

Proposition 3 (INHOMOGENEOUS DIFFERENTIAL EQUATION). For a ma¬ 
trix B G M(n;M) and a vector d G M”' we consider the inhomogeneous linear 
system of ordinary differential equations (with constant coefficients) of sec¬ 
ond order: 

(17) c{t) = Bc{t)-\-d, 

and we consider the one-parameter family of curves Cg defined by Equa¬ 
tion ©• Then we have 

(18) Cs{t) = A{s) c{t)-\-b{s), 
with 

^(s) = ^{1 + cos(s)} 

i.e. Equation Q and Equation ([ 8 |) hold . 

We consider three cases: 

(a) If there is a vector d* such that d = B ■ d^,, then for v,w £ M” the 
unique solution of Equation m with intitial values c( 0 ) = v, c( 0 ) = 
w is given by 

(19) c{t) = c[v, r(;](t) = coB{t) {v + d*) -|- siB(t)(rc) — d* . 
and 

b{s) = ^{cos(s) - 1 } (d*) 

(b) Let B be the nilpotent Jordan block matrix: B = Nn ■= dn(0); cf 
Equation Ii5l Hence Nn{ei) = 0 and Nn{ej) = Cj-i for j >2 for a 
basis d, 62 ,..., Cn o/M”. Then the unique solution c* = c*(t) of 

c;(t) = NnC^{t) + Cn, c*(0) = C*(0) = 0 


j.2n 


( 20 ) 


is given by 
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d G M"' can he written as d = {di, d 2 , ■ ■ ■, dn) = -B • d* + dnCn with 
d* = (0, di, ^ 2 ) ■ ■ ■) dn-i) ■ And the unique solution of Equation (jl7j) 
with intitial values c(0) = v, c(0) = w is given by 

(21) c{t) = coB{t) {v + d*) + siB{t){w) — d* + dnC^{t) 

(22) b{s) = ^ {cob{s) - 1} (d*) + ^ d„c*(s). 

(c) If B = 0 then the unique solution of Equation m with initial values 
V = c{0),w = c(0) and the one-parameter family Cs{t) is given by: 

t^ 1 

c{t) = c[v, w]{t) = d wt V ] Cs{t) = c{t) + - s^d 

and 

b(s) = ^ s^d. 

4 

Proof. The curve 

c{t) = c[v, w](t) = cos(t) (v + d*) + siB(t)(w) — d* 

with c(0) = u,c(0) = w satishes by Proposition [T] 

c(t) = B {coB(t) {v + d*) + siB{t){w)} = B (c(t) + d*) = Bc{t) + B(d*). 

This already proves Equation (I17h in case (a), or rather Equation (I19p . The 
addition rules Equation (|13p show: 

Cs{t) = ^ {c{t — s) + 2c(t) + c{t + s)} 

= ^ [cOBit - s) + COBit + s)] (u + d*) + ^ COB{t) (v + d*) + 

^ [siBit - s) + sis(t + s)] (w) + ^ siB{t){w) - d* 

(23) = i {1 + cos(s)} [coB{t)iv + d*) + siB{t)iw)] - d* 

= Ms){c{t)) +b{s) 

with 

^(•s) = ^ {cob(s) + 1} ; b{s) = ^ {cob(s) - 1} (d*). 

(b) Now let B = Nn- One checks that c^,{t) given by Equation (I20p is the 
unique solution of Equation (fT71) with c(0) = c(0) = 0 and 6 = e^- Then for 
any solution c = c{t) of Equation (|17p the curve c{t) — dnC^.{t) is a solution 
of Equation 

e{t) = B (e(t) + d*) = Be{t) + (di,..., d„_i,0) . 

Hence we conclude from Part (a) that 

e{t) = COBit) (v + d*) + siB{t)iw) - d*. 
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This implies that c{t) satisfies Equation (I2ip . Since 

I {c*(t - s) + 2c*(t) + c*(t + s)} = i {cob(s) + 1} (c*(t)) 

together with Equation (f23l) proves Equation Q and Equation (|22p . 

(c) For B = 0 Equation (|17p is c{t) = d, hence we obtain c{t) = dt'^/2+wt+v . 
Furthermore cob{s) = 1 = A(s) and Cs{t) = ^ {c{t — s) + 2c{t) + c{t + s)} = 
c{t) + \s‘^d. □ 


Remark 3 (Wave equation). Let c = c{t) be a smooth curve with the 
one-parameter family of curves Cg defined by Equation ([5]). Then 


= 1 eft - 
ds'^ 4 

Then Cs{t) := Cs{t) - ^ c(t) 
the wave equation: 


with initial conditions Co(t) 


1 d'^Cgit) 

s) + I c(t + s) = 

= {c{t - s) -I- c{t + s)} /4 


d'^Csjt) _ d‘^Cs{t) 

ds^ dt"^ 

= c(t)/2 and 


dehnes a solution of 


dcs{t) 


If the smooth curve c is a solution of Equation (I17p then Proposition [3] 
implies: 


Cs{t) = (^(s) - 1/2) (co(t)) + h{s) = cob{s) (co(t)) + b(s). 

Hence during the evolution s i—)■ c<j the afhne form of the curve c = cq is 
preserved. This motivates the notion soliton for these solutions of the wave 
equation. 


In the sequel we study which invertible matrices D can be written in the 
form (1 + COB{s)) l‘2‘ for some real matrix B. Since for any invertible matrix 
U we have cojjBU-^i^) — Ucob{s)U~^ it is sufficient to check the possible 
Jordan normal forms of cob{s) for the different Jordan blocks B. It turns 
out that a large class of invertible matrices D can be written in this form. 

Proposition 4 (The image of {B,s) i—)• (1-|- cob(s)) /2). Consider the 
mapping f : {B,s) G M(re;M) x M i-A f{B,s) := (1 +cob(s)) /2. An invert¬ 
ible matrix D G M(n;]R) lies in the image of f if and only if the (real or 
complex) Jordan normal form of D satisfies the following: If X < 0 is a real 
and negative eigenvalue of D and if for some m > 1 the Jordan block JmW 
occurs in the Jordan normal form, then the number of Jordan blocks JmW 
in the Jordan decomposition of D is even. 

In particular a diagonizable matrix D can be written as: D = f{B,s) = 
(1 + cob{s)) /2 for some [B, s) G M(n;M) x M if and only if the eigenspaces 
of all real and negative eigenvalues are even-dimensional. 
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Proof. We compute for the possible complex Jordan normal forms J of a real 
matrix B the complex Jordan normal form of coj(s) resp. /(J, s). Since the 
matrix D = (1 + cob('S)) /2 is supposed to be invertible we exclude —1 as 
eigenvalue of cob (s). 

(a) If J = A • 1, A € M, then f{J,s) = (1 + cosa(s))/2 • 1. Note that 

{(1 + cosa(s))/2 ; A, s G M} = {t G > 0}. Hence f{J,s) is a diago¬ 

nal matrix with non-negative real eigenvalues. 

(b) Assume the matrix B contains a complex Jordan block J = u ■ Im; u G 

C \ M, then also J = u • 1 is a (distinct) Jordan block, i.e. there is an even¬ 
dimensional invariant subspace on which the Jordan normal form is given 
byj0j = u ■ Im ■ Im- Then / (J © J, s) = (1 + cosh(rcs)) /2 • 1 © 
(1 + cosh(u;s)) /2 • 1, here w £ C,u = . Note that cosh : C —^ C is surjec¬ 

tive. If cosh(tcs) G M the corresponding eigenspace is even-dimensional. 

We conclude from (a) and (b): A real and invertible matrix D which is 
diagonalizable over C can be written in the form D = f{B, s) for a matrix 
B diagonalizable over C and a real number s if and only if it does not have 
a real and negative eigenvalue A with an odd-dimensional eigenspace. 

(c) Now assume that J = Jm{u) with m > 2. For u = 0 ; we obtain for 
Nm = Tm(0) : 

2 2m—2 

(24) coiv^(s) = Im + + • • • + ^ ■ 

Note that the matrix satisfies N^{eiJ^k) = ep / = 1,2,..., m—k, N^{ei) = 
0, Z = 1, 2,..., fc — 1. Since for s / 0 


rank (coAT^ (s) - 1 ^) = m - 1 


we obtain: For s 7 ^ 0 the matrix con^{s) is conjugate to the Jordan block 
Jm(l) and f(Nm,s) is also conjugate to the Jordan block Jm(l). 

For u 7 ^ 0 we obtain: 


00 2k 

(25) = E ' ■ '"‘-'iVi 


- ^2k 


k=0 


(2k)l 

k=0 ^ ’ j=0 


3 


2k 

* k 


E 

k=0 


{2k)\ 


• Im + X^ 


„ 00 „2fc-l 

* ^ ^ h — A 

2 ^ ^ + 


+a2[s)N‘^ + ... + am-i{s) ■ N, 


m—1 

m 


for some 0 : 2 ( 5 ), ■■■, Om-i(s) G C. If u = A G M we have 


coj(s) = cosa(s) • Im + I sinA(s) • Nm + a2N^ + ... + Um-lNff ^ . 

Note that {cosa(s) ; A,s G M} = {x G M, x > —1} and that sinA(s) 7 ^ 0 
whenever cosa( 5 ) 7 ^ ±1. Hence for cosa( 5 ) 7 ^ ±1 the matrix coj(s) satisfies: 


rank (coj(s) — cosa(s) • Im) = m, — 1 
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i.e. the matrix / (Jm(A), s) is conjugate to the Jordan block ((1 + cosa(s)) /2) . 
Thus we have shown that a Jordan block with real fJ. ^ 0 and m > 2 

can be written in the form / ((A), s) for some real A and some s G M up 
to conjugacy if and only if /r is positive. 

(d) Now assume that the real matrix B contains a complex Jordan block 
J = Jmiu),u G C \ M. Then let re G C \ M be a square root, i.e. u = uP'. 
Hence also J = Jm{u) is a (distinct) Jordan block of B, or rather there is 
an invariant subspace on which the Jordan normal form of A is given by 
J (B J = Jm{u) (B Jm(u). We conclude from Equation [25l 

coj(s) = cosh('u;s)' ^ smh{ws)Nm + a 2 {s)N^ + ... + am-i(s) • 

for some a 2 ('S), • • •, am-i{s) G C. If cosh(t(;s) / ±1 we have sinh(t(;s) p 0, 

i.e. 

rank (coj(s) — cosh(r(;s) • Im) = m — 1 

or rather the matrix coj(s) is conjugate to a Jordan block Jm (cosh(i(;s)) and 
/ (j © j) is conjugate to Jm ((1 + cosh(t(;s)) /2) © Jm ((1 + cosh(u;s)) /2) . 

We conclude from (c) and (d) that an invertible and real matrix D whose 
Jordan normal form contains a Jordan block J = Jm{v),m > 2 can be 
written in the form f{B,s) for a real matrix B and s G ffi if and only if 
the following holds: If v is real and negative u = A < 0 then the number of 
Jordan blocks Jm (A) in the Jordan decomposition of B is even. □ 

3. Discrete Curve Shortening 

An (infinite) polygon x = in MA is defined by its vertices Xj G M”. 

We call V = P(M”) the vector space of these polygons. We can identify the 
polygon X with the piecewise linear curve x : M —)■ which is a straight 

line on any interval [j,j + 1] and satisfies x{j + u) = (1 — u)xj + uxjj^i for 
any u G [0,1], j G Z. If there is a positive number N such that xj+at = Xj 
for all j G Z, then we call the polygon x closed or rather periodic with N 
vertices or rather of period N. In this case we can identify the index set with 
Zjv = Z/(A^ • Z). We denote the set of closed polygons with N vertices by 
Vn = The midpoint mapping is given by 

M : X G P(M”) ^ M(x) G iP(M"); (M(x)),. := \ (xj + xj+i). 

For a closed polygon x G Pjv(M”) its length is given by 
T (t) := ~ ®ill > here ||.|| denotes the Euclidean norm. 

The triangle inequality implies the following curve shortening property of 
the midpoint mapping in the general case: 

m 

L{M{x)\[j,j + m]) = J]] ||(M(x))j.+;. - (M(x))j.+;i,_i 

k=l 

< L{x\[j + 1/2, j + m + l/2]) 
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with equality if and only if the points Xj, Xj+i,..., Xj+m+i lie on a 

straight line. Thus the midpoint mapping M is length decreasing on closed 
polygons, i.e. for all closed polygons x we have: L(M(x)) < L(x) with 
equality if and only if x is constant, i.e. L(x) = 0. We modify the midpoint 
mapping as follows: 

Definition 2 (Curve shortening process). We introduce the following 
mapping T : P(M”) : 

{T{x))^ = j {xj-i + 2xj + Xj+i} . 

(a) We call a polygon x = {xj)j£z affinely invariant under the mapping T 
or rather an (affine) soliton of the curve shortening process T) if there is an 
affine map {A, b), A £ Gl , b € M” such that for all j £ Z : 

(26) T{x)j = ^ {xj-i + 2xj + Xj+i} = A (xj) + b. 

(h) We call a smooth curve c : M ^ M” affinely invariant under the mapping 
T (or rather an (affine) soliton of the curve shortening process T) if there 
is some e > 0 such that there is a one-parameter family s £ (—e, e) i-A 
{A{s),b{s)) £ G/(n,M) x M"' of affine maps such that 

(27) Cs{t) = ^ {c(t - s) + 2c{t) + c{t + s)} = A{s)c(t) + b{s) 
for all s G (0, e),t G M. 

It is obvious that these notions are affinely invariant. For a G M, s G (0, e) 
the polygon x = x{a, s) with x(a, s)j = c{a + sj),j £ Z lying on a smooth 
curve c : M —M” is a soliton of the curve shortening process T if the curve 
c is also a soliton of the corresponding process T on curves. On the other 
hand: A polygon x = which is a soliton of the curve shortening 

process T satisfying Equation (1261) can be obtained from a smooth curve, 
which is a soliton as defined in Equation (|27p if and only if A can be written 
in the form A = (1 +cob(s)) /2. In Proposition 0] the Jordan normal form 
of these matrices are classified. 

Remark 4 (Eigenpolygons of T). If we consider closed polygons then 
the midpoint mapping defines a linear map M : Vn —^ 'Pn on the (n • 
A^)-dimensional vector space Vn, and one can use a decomposition into 
eigenspaces, cf. HU and [2]. The matrix is in particular circulant. Eor 
n = 2 one can identify the complex numbers C with M^, then the eigenval¬ 
ues are = (1 -|- exp(27riA:/A^)) /2 , A: = 0,1,... , — 1. The corresponding 

eigenvectors are the polygons 

(28) ..., G ; zj^^ = exp {2'Kijk/N) 
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for k = 0,1,..., N — 1. Then the linear map T : Vn —^ 'Pn is also circulant 
and has the form 


/ 2 100... 000l\ 

1 2 1 0 ... 0 0 0 0 

0 1 2 1 ... 0 0 0 0 


T = 


1 

4 


0 0 0 0 

0 0 0 0 

\ 1 0 0 0 

with eigenvalues 


1 

fJ-k = fJ-N-k = 2 



12 10 
0 12 1 
0 0 12 / 


= 0,1,... ,1V- 1 


and eigenvectors given by Equation (1281) . Note that all polygons 
given by Equation (|28l) lie on the unit circle c{t) = exp(27rit) and = 
c{kj/N) . The eigenpolygons are solitons on which the map T is a ho- 
mothety. Although the curve c is simple and convex the polygons form 
regular A^-gons which are simple and convex only for A; = 1, — 1, cf. [21 

Fig. 5]. 


Proposition 5 (assignment matrix A polygon). Let (A,b) : x £ 
M” I—>■ Ax + b G M"' be an affine map and u,v £ be two points in M”, and 
jo £ Z. Then there is a unique polygon x £ 'P{MT) with Xj^^ = u, = v 

which is affinely invariant (with respect to A and b) under the mapping T. 


Proof. If X € P(M”') is affinely invariant under the mapping T we have 
(29) {T{x))j = ^ {xj-i + 2xj + Xj+i} = A (xj) + b. 

Hence the sequence (xj) with xj^ = u, = u is uniquely determined by 

the recursion formulae 


(30) Xj_|_i = 2 {2A — 1) (xj) —Xj_i+46; Xj_i = 2 {2A — 1) (xj) —Xj+i+46. 

□ 


For a given smooth curve c : M —M” we define the one-parameter family 
Cs : M —)• M"' by Equation ([5|). The curves Cg are obtained from c = cq by 
applying the mapping T as follows. For arbitrary o G M, s > 0 let x = x(a, s) 
be the polygon Xj = Xj{a, s) := c (a + js) , j G Z. Then 

Cs (o + js) = {T{x))j = i {c(a (j - l)s) + 2c(o -h js) + c{a + (j + l)s)} . 

Hence the vertices (T (x(a, s)))j of the image T (x(a, s)) of the polygon x = 
x{a,s) under the mapping T lie on the curve Cg, or rather the curve Cg is 
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formed by the images T [x (a, s)) of polygons of the form x = x(a, s) on the 
curve c. 

Theorem 1 (SOLITONS AS SOLUTIONS OF AN ode). Let c : M — s- M"' be 
a smooth curve such that the one-parameter family defined by Equation ([5]) 
defined for all s £ (—e, e) and some e > 0 satisfies: 

Cs{t) = A{s){c{t)) + b{s) 

for all t £ M,—€ < s < e and for a smooth one-parameter family s i-A 
A(s) £ Gl{n,M) of linear isomorphisms and a smooth curve s e-)• b{s) £ M”. 
I.e. the curve c is affinely invariant under the mapping T, cf Definition\^ (b). 
Assume in addition that for some tg £ K the vectors c{to),c{to ),..., (tg) 
are linearly independent. 

Then the curve c is the unique solution of the differential equation 

(31) c{t) = Bc{t) + d 

with initial conditions v = c{0),w = c(0). The functions s e-)• ^(s) £ 
M (n,M), s I—)• b{s) £ M"" satisfy the differential equations 

(32) A'fs) = {A{s) - 1/2) 5; b''{s) = {A{s) - 1/2) d 

with B = 2A"(0), d = 26"(0) and initial conditions ^(0) = 1, 6(0) = 0, ^^(0) = 
0, 6'(0) = 0. In particular, we have A{s) = (1 + cob{s)) /2, cf. Proposition\^ 
The possible Jordan normal forms of the matrices ^(s) are given in Propo¬ 
sition^ 

For a one-parameter family s e-c* of curves or a one-parameter family 
s e-)• ^(s), s e-)• 6(s) of matrices ^(s) or vectors 6(s) we denote differentiation 
with respect to s by Differentiation with respect to the curve parameter t 
of a one-parameter family of curves Cs{t) or a single curve t i-A c{t) is denoted 
by c. 

Note that for the potential 

(33) U{x) = —-{Bx,x) — {d,x) 

we can write instead of Equation (fSTI) the following Equation: 

(34) c{t) = —grad 6/ (c(t)) . 

It follows that the function 

fit) ■= l|c(t)f - ^{Bc{t),c{t)) - {d,c{t)) 
is constant for any solution of Equation (1311) . 

Proof. Let 

(35) Cs{t) = A{s)c{t) -\- b{s) = ^ {c{t - s) + 2c(t) -|- c{t -\- s)} . 
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For s = 0 we obtain c{t) = Co{t) = ^(O)c(t) + b{0) for all t G M, or rather 
(A(0) — 1) (c(t)) = —6(0) for all t. We conclude that 


(36) 


(^(0) -1) (c(^)(t)) =0 


for all k > 1. Since for some to the vectors c{to), c{tQ),..., (to) are lin¬ 
early independent by assumption we conclude from Equation (I36p : ^(0) = 
1, 6(0) = 0. Equation (l3^ implies for A; > 1 : 


and hence 


A{s)c^'^\t) = i {c(^)(t - s) + 2c^^\t) + + s)} . 

A\s)c^'^\t) = |c(^+i)(t - s) - + s)} 


or rather 

^'(0)c(^)(t) = 0. 

Since for some to the vectors c(to)) ■ ■ ■ are linearly independent we 

obtain A^(0) = 0 and hence by Equation (l3^ : 6'(0) = 0. 

We conclude from Equation (f35ll : 

d'^Csit) 


ds"^ 


= A"{s)c{t) + h"{s) 


Since ^(0) = 1 the endomorphisms A(s) — 1/2 are isomorphisms for all 
s G (0, e) for a sufficiently small e > 0. Hence we obtain for s G (0, e) : 

(37) c{t) = (A(s) - l/2)-i A"is)c{t) + (H(s) - 1/2)-^ b"{s). 

Differentiating with respect to s : 


(8)-l/2)-'2l"(»))'(c(t)) + 
and differentiating with respect to t : 

is)-l/2)-^A"{s)yc^^\t) = 0-,k = 0,l,2,..., 


n . 


Since c{to), c{to), ■ ■ ■, (to) are linearly independent by assumption we con¬ 
clude that 

((21(s)-l/2r'A"(s))' = 0. 

Hence with B = 2H"(0),d = 26"(0) we obtain 

H"(.) = (H(s) - 1/2) B ; b"{s) = {A{s) - 1/2) (d), 
i.e. the differential Equations (l3^ . Eor s = 0 we obtain from Equation (l371) : 

c(t) = Bc{t) + d , 

i.e. Equation (l3T]) . □ 
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We can combine the results of Theorem [T] and Proposition [3] to give the fol¬ 
lowing characterization of affinely invariant curves under the affine mapping 
T as solutions of an inhomogeneous linear differential equation of second 
order: 

Theorem 2 (characterization of affinely invariant curves as 
SOLUTIONS OF AN ODe). (a) Let c : M —)• M” be a smooth curve affinely 
invariant under the mapping T. Assume in addition that for some to € M 
the vectors c(to)) c(to)) • • • > linearly independent. Then there is a 

unique matrix B S M(n;M) and a unique vector d £ M” such that 

c{t) = Bc{t) + d. 

(h) Let B be a real matrix and d a vector in M”. Then any solution c = c{t) of 
the inhomogenous linear differential equation c{t) = Bc{t) + d with constant 
coefficients defines an affinely invariant smooth curve under the mapping T. 

For a closed polygon x G Vn the center of mass Xcm is given by 

1 ^ 

(38) Xcm = Xj ■ 

i=i 

Since (T(x))^^ = Xcm we conclude: There is no translation-invariant closed 
polygon, since the center of mass is preserved under the curve shortening 
processes B or rather T. 

Remark 5 (Generalization of the map T). For three points x,y,z G 
M"' define the affine map 

r : X X M" ^ ; T{x, y,z) = ^{x + 2y + z}. 

Hence the mapping T : V{MT) —)• P(M"’) introduced in Definition [2] satisfies 
for all j G Z : 

T{x)j = T {xj-i,Xj,Xj+i) . 

The one-parameter family Cg associated to the smooth curve c by Equa¬ 
tion Q can be written as: 

Cs{t) := T (c(t - s),c{t),c{t -I- s)) = ^ {c(f - s) -I- 2c{t) + c{t + s)} . 
fn the following we will allow a slightly more general curve shortening process 

(39) r„ : iP(M") ^ V{Rn, {Ta{x))j := (x,_i, x,-, x,+i) 

based on the affine map Ta ■ ^ 1^” ]Ta{x, y, z) = ax + (1 — 2 a)y + az 

for a ^ 0, i.e. T = T 1 / 4 . For a = 1/3 the point Tii 3 {x,y, z) is the center 
of mass {x + y + z)/3. The curve a G [0,1/2] i-)- T/* (xj_i, Xj, Xj+i) G M"" 
is a parametrization of the straight line connecting Xj with the midpoint 
(xj_i +Xj^i)/2 of the points Xj-i,Xj^i. These mappings Tq, are considered 
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for example in [21 p.238-39] and |3l ch.5.1]. For a smooth curve c one defines 
the associated one-parameter family of curves 

Ca,s{t) = Ta {x{t — s),c{t),c{t + s)) = ac{t — s) -|- (1 — 2 q:) c{t) + ac{t + s). 

We call a smooth curve c affinely invariant (or a soliton) with respect to if 
there is a one-parameter family (^q,(s), &a(s)), s £ (—e, e) for some e > 0 of 
affine mappings such that Caffit) = Aa{s)c{t) + ba{s) for all t £ M, s £ (—e, e). 
Then 

Cs{t) = Ci/4^sit) = c{t) + ^ {c{t - s) - 2c{t) + c{t + s)} 

_ (4a - l)c(t) Cq^sjt) 

4a 

and 

Ca,s{t) = 4aCs(t) - (4a - l)c(t). 

We conclude: A smooth curve c is a soliton for the transformation for 
a 7 ^ 0 if and only if it is a soliton for the transformation T = T ]^/4 and for 
the corresponding affine maps {Aa(s),ba{s)) we obtain: Aa(s) = 4aA(s) — 
(4a — 1)1 ; ba{s) = 4a5(s). 

4. Semidiscrete flows of polygons 

The mapping T introduced in Definition [2] or Tq, defined in Remark [5] or 
Equation (1391) can be seen as a discrete version of the semidiscrete flow 
defined on the space P(M”) of polygons: For a given polygon the 

flow s e->■ {xj{s)) £ P(M”') is defined by the Equation 

(40) ^ - ‘2.Xj{s) + Xj+i(s) , Xj(0) = Xj . 

This fiow is discussed for example in |5]. ft is a linear first order system of 
differential equations with constant coefficients forming a circulant matrix. 
Hence one can write down the solutions explicitely. If we approximate the 
left hand side of this equation by {xj{s -|- a) — Xj{s)) /a we obtain 

Xj{s + ot) = axj-i{s) + (1 — 2a) Xj{s) + axj+i{s) 

= Ta{Xj-l{s),Xj{s),Xj + l{s)) . 

Therefore the mappings T and Ta can be seen as discrete versions of the 
flow Equation (jdOp . fn |5] the flow s e-?> Xj{s) introduced in Equation (I40p 
is called semidiscrete since it is a smooth flow on a space of discrete objects 
(polygons). Then a discretization of the semidiscrete flow yields to the dis¬ 
crete process T and Ta discussed here. On the other hand the connection 
of the semidiscrete flow with the smooth curve shortening flow in Euclidean 
space is discussed in detail in |5l Sec. 5]. 
ff we consider the functional 
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on the space Vn of closed polygons then we obtain for a curve s G (—e, e) e-)• 
x{s) = with X = x( 0 ) and x = i;( 0 ) : 

N 

= - ^ {xj,Xj-i - 2xj + Xj+i) 
j=0 

and we obtain for the gradient grad-F 2 ( 3 ^) : 

gradF 2 (x) = - (xj-i - 2xj + Xj+i)j^^^ . 

Hence the semidiscrete flow can be viewed as the negative gradient flow of 
the functional F 2 , cf. [5l Sec. 6 ]. An affine transformation x G M” e-)- ^(x) = 
A{x) + 6 G MA on MA induces an affine transformation A on Vn ■ A = 
(xi,... ,Xn) = (xi),...,^ (xat)) . In contrast to the functional F 2 its 
gradient gradF 2 is invariant under A : 

gradF 2 = - {A (xj_i) - 2A {xj) + A (xj+i)).^^^ = A{gradF 2 (x)) . 

Definition 3 (affine invariance under the semidiscrete flow). We 
call a smooth curve c : M —)• M"' affinely invariant under the semidiscrete flow 
given by Equation iO]) (or a soliton ) if there is some e > 0 such that for any 
s G (0, e) there is an affine map such that the one-parameter 

family 

(41) Cs{t) = A{s)c{t) + 6 (s) 
is a solution of the flow equation 

= Cs{t - 1 ) - 2Cs{t) + Cs{t + 1 ) 

os 

for a// s G (0, e), f G M. 

In the following Proposition we show that the solitons of the mapping T or 
Tq coincide with the solitons of the semidiscrete flow given by Equation (j40]) : 

Proposition 6 (affinely invariant curves under the semidiscrete 
flow). Let B be a matrix and d a vector in M” and let c = c{t) be a 
solution of the inhomogenous linear differential equation c(f) = Bcif) + d 
with constant coefficients for which c( 0 ), c( 0 ),..., ( 0 ) are linearly inde¬ 

pendent. Then the curve c defines an affinely invariant smooth curve under 
the semidiscrete flow given in Equation dini) . 

Proof. We conclude from Proposition [3] and Equation (fT 8 )l that there is a 
matrix Ai = (A(l) — l)/4 and a vector bi = b(l)/4 such that 

(42) c{t — 1) — 2c{t) + c(t + 1) = Aic{t) + 61 
for all f G M. Let 

A{s) = exp{Ai ■ s) ] b{s) = / exp(Ai • cr) ( 61 ) dcr, 

Jo 


dF2{x{s)) 

ds 
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Figure 3. Scaling, 
Case (la) 



Figure 4. Scaling, 
Case (lb) 


and 

(43) Cs{t) = A{s)c{t) + b{s). 

Using Ai ■ ^(s) = A{s) ■ Ai we obtain: 

= A'{s)c{t) +V{s) = Ai • A{s)c{t) + Ais){h) 

= A{s) {Aic{t) + bi) . 

We conclude from Equation (|42|) and Equation ([43]): 

= 2(s) (c(t - 1) - 2c{t) + c{t + 1)) = Cs{t - 1) - 2Cs{t) + Cs{t + 1). 
os 

Hence the curve c is affinely invariant under the semidiscrete flow. □ 


5. Planar solitons 

We study the planar case n = 2. We conclude from Theorem [2] that solitons 
c are solutions c{t) = {x{t),y{t)) of the differential equation 

c{t) = Bc{t) + d. 

We discuss these solutions using Proposition [2] and Proposition [3j 

Let B S M{2; M) be a matrix in (real) Jordan normal form. Then we consider 

the following cases: 

(1) B is diagonalizable (over M) and invertible, and d = 0. i.e. 

5= ( 0 ) , 61,62 GM-{0}, 

this affine map is called scaling. If the diagonal entries coincide the trans¬ 
formation is also called a homothety. The differential equation c = Be 











20 


CHRISTINE RADEMACHER AND HANS-BERT RADEMACHER 



Figure 5. Spira 
mirabilis, Rotate 
and scale, Case ( 2 a) 



Figure 6 . Rotate 
and scale, Case (2b) 


implies x = bix, y = b 2 y- With the notation introduced in Equation ([9]) 
we obtain with c( 0 ) = (i;i, ^ 2 ), c( 0 ) = {wi,W 2 ) : 


c{t) = {vi cosftj {t) + wi siuftj {t),V 2 COS 62 {t) + W 2 sinb 2 (t)) 


The matrices ^(s) are also diagonal matrices (or scalings) 




1 


{1 + cob{s)} 


1 

2 


1 + COSb^(s) 
0 


0 ^ 

1 + COSb 2 (s) J 


In particular the diagonal entries of the matrices ^(s) are non-negative 
for all s. Hence the one parameter family Cg of curves is produced by 
scaling (with a diagonal matrix) from the soliton c = cq. 

Particular examples are: 

(la) We obtain closed curves if 61 = —< 0, 62 = — < 0 and A 1 /A 2 G Q. 

Then we obtain for example Lissajous curves of the form 


c{t) = (rci sin (Ait), V2 cos (A 2 t)) 

see Figure [ 3 ] for Ai = 4, A 2 = 9,r(;i = V 2 = 1 ,W 2 = vi = 0 and 
0 < t < 6.3. Another example is given in the Introduction, it is the 
curve c(t) = (cos (2t), cos (3t)) . In Figure[2]we show the curve c = c{t) 
and C 0.4 = co. 4 (t) = A (0.4) c{t) : 

coA{t) = 2 (0-8)) cos(2t), (1 -|- cos (1.2)) cos(3t)) . 

(lb) If 61 = A^ > 0,62 = —^2 < 0 we obtain for example curves of the 
following form: 


c(t) = (rci sin (Ait), V 2 cosh (A 2 t)) 

In Figure[4]this curve is shown for the parameters Ai = 8 , A 2 = 1, V 2 = 
wi = l,vi = W 2 = 0 and —1.3 < t < 1.3. 
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(Ic) If B is diagonalizable and has the form 


A = 


0 0 
0 62 


,^2/0. 


and d = 0 then x(t') = 0 ,y = b 2 y, i.e. c{t) = (t, V 2 COS 62 (t) + W 2 sinb 2 (t)) 
are solutions, for example 


c{t) = (t, exp(f)), c(t) = (t, cosh(t)), c(f) = (f,sin(t)). 

(2) The following case corresponds to the case that the matrix B has no real 
eigenvalues. Hence H is a similarity, i.e. a composition of a rotation and 
a dilation x Xx for some A 7^ 0 . We identify with the complex 
numbers C and assume that the matrix A is complex linear, i.e. can 
be identified with the multiplication with a non-zero complex number y, 
Then we are looking for a solution 2 : : t G M 1-^ z{t) = x{t) + iy{t) G C 
of the differential equation z = yz. For a complex number w with y = 
up', lu = ui + iu2,ui,U2 G M,= 'R{'w),U2 = ^{w) a solution has the 
form 

z{t) = hi exp{wt) + /i 2 exp{—wt). 

for /ii,/i2 G C. If we write hi = hu + f/112,/12 = /121 + 2/122 with 
hii,hi2,h2i,h22 G M then we obtain for x{t) = ?fi{z{t)),y{t) = ^{z{t)) : 

x{t) = {/ill exp(uif) +/i2i exp(—uif)} cos(u2/) 

+ {-hi2 exp(uit) + /122 exp(-uif)} sin(u2/) 

y{t) = {/112 exp(uit) +/122 exp(-uit)} cos(u2/) 

( 44 ) + {/ill exp(riit) — /121 exp(—riit)} sin(M2/) 

For arbitrary hi , /12 the one-parameter family Zg (t) in complex notation 
is given by 

Zs{t) = ^{z{t - s)+ 2z{t) + z{t + s)} 

= 2 (^1 ex.p{wt) + /i2 exp(—rut)) . 

In real notation we obtain the matrix: 

.. . _ 1 / 1-|-cosh(uis) cos(u2s) — sinh(uis) sin(M2s) 

^ 2 \ sinh(uis) sin(tt2s) 1-|-cosh(iiis) cos(u2s). 


Particular examples are: 

(2a) The logarithmic spiral [spira mirabilis): 

c{t) = exp(tiit) (cos(ti 2 /), sin(u 2 t)) . 

i.e. /ill = 1)/2 i2 = h 2 i = /122 = 0. Figure [5] shows this curve for 
ui = 0.3, U 2 = 4 and —3 < / < 3. 


22 


CHRISTINE RADEMACHER AND HANS-BERT RADEMACHER 



Figure 7. Rotate 
and scale, Case (2c) 


Figure 8. Shear, 
Case (3) 


(2b) The curve 


c{t) = (cosh(riit) cos(ri2t), sinh(Mit) sin(rt2t)) 


(i.e. /ill = /i2i = 1/2,/ii2 = /i22 = 0) is shown in Figure [6] for the 
values Ml = 1, M2 = 20 and 0 < / < 1.2. 

(2c) The curve given in Equation (1441) for hn = 1,/121 = 1.3,/112 = /122 = 
0,Mi = 1,M2 = 20 is shown in Figure [7] for —0.57 <t< 0.885. 

(3) Let B be non-zero and nilpotent, i.e. 



and d = { 0 ,d 2 )- 

Then a solution is given by c(/) = (c/2/^/24-|-os/^/G + a2t^/2-|-ai/-|- oq,/) 
for aQ,ai,a 2 ,a 3 £ M. The matrices ^(s) are of the form: 



cf. the Proof of Proposition 01 Hence the one-parameter family s i-A 
is formed by shear transformations. The curve with parameters c/2 = 
0.1,03 = 0.2,02 = —4, oi = —1,00 = 0 for —30 < / < 25 is shown in 
Figure El 

(4) Let B be invertible with real eigenvalue and not diagonizable, i.e. 



with b £ W,b 0 and d = 0. Then c{t) = {x{t),y{t)) with 



y{t) = V 2 COSb{t) + W 2 smb{t) 
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Figure 9. Shear 
and Scaling, Case (4) 



Figure 

10. Translation, 
Case (5) 



Figure 11. Translation and Scaling, Case (6) 

The matrices A(s) are of the form 

,.1/ l + cosb(s) s-smb{s)/2\ 

~ 2 V 0 1 + cosb(s) J ’ 

cf. the Proof of Proposition 01 Hence the one-parameter family s i—)■ c* is 
formed by a composition of shear transformations and scalings. For the 
parameters b = —1,111 = l,i;2 = —0.1,t(;i = —10, r(;2 = 1 and —30 <t< 
40 the curve is shown in Figure [9l 

(5) If H = 0 and d ^ 0, then we obtain (up to an affine transformations) the 
parabola c{t) = {t,t^) as translation-invariant curve, cf. Proposition [3Kc) 
and Figure [TOl 

(6) If B is of the form 


B = 


0 0 
0 b 
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with non-zero b and b = (2,0). Then xi{t) = 2 and X 2 {t) = bx 2 {t). Then 
c{t) = cosb{t) + wi sinfe(t)) 

Examples are (for b = ±A^ :) 

c{t) = (f^,cos(At)) ,c{t) = (t^,exp(At)) ,c(t) = cosh(At)) . 

In Figure nn the curve c(t) = (t^,sint), —10 < t < 10 is shown. 
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